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O . Abstract 
O 

^ We present predictions for flavour- violating charged-lepton decays induced by the seesaw 
^ mechanism implemented within the constrained minimal supersymmetric standard model 
^ . (CMSSM) with universal input soft supersymmetry breaking terms. We assume that one 
j> [ heavy singlet neutrino almost decouples from the see-saw mechanism, as suggested by 
I the pattern of light neutrino masses and mixing angles. This is suggested independently 
^ • by sneutrino inflation with a low reheating temperature, Thh ~ 10^ GeV, so as to avoid 
■ ^ > overproducing gravitinos. This requirement further flxes the mass of the weakly-coupled 
sneutrino, whose decays may lead to leptogenesis. We flnd that BR^jj, e'j) > 10^^^ 
but BR{t fi'y) < 10^^ in the bulk of the acceptable parameter space, apart from a 
few isolated points. The ratio BR{n e'y)/BR{T 67) depends on only one complex 
parameter, and is particularly interesting to compare with experiment. 



1 Introduction 



The observation of neutrino mixing |Tj has led many authors to consider the possibility 
of flavour violation and CP violation in the charged-lepton sector [2]. The most predic- 
tive framework for such studies is the seesaw model for neutrino masses ^ implemented 
within the minimal supersymmetric extension of the Standard Model, constrained to have 
universal soft supersymmetry breaking terms at some input renormalization scale charac- 
teristic of grand unification (CMSSM). This framework allows one to explore the possible 
links between lepton-flavour violation (LFV), neutrino oscillations and leptogenesis [4^ 5j, 
assuming the cosmological baryon asymmetry to have originated from CP-violating decays 
of heavy singlet neutrinos and their supersymmetric partners. 

Even the minimal three-generation seesaw model contains 18 free parameters 0111, in 
addition to those describing supersymmetry breaking. Four combinations of these seesaw 
parameters have been determined by neutrino-oscillation experiments, two mixing angles 
and two mass-squared differences. More of the seesaw parameters may be determined 
by low-energy neutrino experiments, but not all of them, and speciflcally not all those 
controlling leptogenesis |H1 13 UHl HH Il2|- Under these circumstances, supplementary 
hypotheses are needed if one is to make unambiguous predictions for LFV, leptonic CP 
violation and leptogenesis. Fortunately, one can find good physical motivation for certain 
simplifying assumptions. The observed hierarchy of differences in neutrino masses-squared 
can be interpreted as a hierarchy in the masses themselves. Such a pattern of the active 
neutrino masses together with the measured nearly bi-maximal neutrino mixing would 
be naturally explained if one of the heavy singlet neutrinos is almost decoupled from 
the seesaw mechanism [121 El . As we analyze in this paper, this decoupling hypothesis 
imposes important constraints on the seesaw parameters, and hence leads to interesting 
predictions for LFV processes. 

The decoupling hypothesis is supported by the suggestion that the scalar field supposed 
to be responsible for infiation, the infiaton, is one of the heavy singlet sneutrinos fTS]. 
Requiring a low reheating temperature, Thh ~ 10^ GeV after infiation, in order to solve 
the cosmological gravitino problem [1^, forces the infiaton sneutrino to couple very weakly 
to ordinary matter and its superpartner almost to decouple from the seesaw mechanism. In 
order to explain simultaneously the duration of the infiationary epoch responsible for the 
large-scale structures observed in the Universe and the magnitudes of the perturbations 
observed in the cosmic microwave background, the mass of this sneutrino should be around 
2 X 10^^ GeV This is well inside the range thought plausible in the seesaw model. 
Also, as discussed elsewhere, sneutrino infiation makes predictions for cosmic microwave 
background (CMB) observables, such as the scalar spectral index and the ratio of tensor 
to scalar perturbations pTJ, which are consistent with data on the CMB and cosmological 
structure formation 

It is striking that, on the one hand, naturalness arguments for models of neutrino 
masses and mixings and, on the other hand, the hypothesis of the sneutrino infiation, 
both motivate independently a similar pattern in the seesaw mechanism. 
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The scenario with the heaviest right-chiral neutrino decoupled is a natural possibility. 
However, a hierarchy of the heavy singlet neutrino masses could be compensated by 
some hierarchy in the neutrino Yukawa couplings, and it is interesting to consider also 
the decoupling of one of the lighter neutrinos. In fact, decoupling of the lightest singlet 
neutrino was considered in the sneutrino inflation model of [Ej, but one may also consider 
inflation driven by one of the heavier sneutrinos, if it is decoupled. The purpose of this 
paper is to outline the different scenarios for decoupling one of the heavy singlet neutrinos, 
and explore in some detail the LFV signature of some specific examples. 

We discuss our different decoupling assumptions for the neutrino Yukawa couplings in 
Section O and explore their various predictions for LFV decays in Section [HI For each case 
we also discuss the possible leptogenesis, which is necessarily nonthermal in the sneutrino 
infiation models with low reheating temperature, with the CP asymmetry generated in 
the decays of the infiaton |Tn|. As we discuss in Section [3 the details of this mechanism 
depend on which neutrino plays the role of the infiaton. 

The main result of our investigations is that the decoupling hypothesis leads to rela- 
tively rigid predictions for the branching ratios BR{ii —>■ ij'j). For most of the acceptable 
parameter range, we predict BR{fi 67) > 10~^^, within the reach of the experiment 
now underway at PSI. However, for generic models with BR{fi 67) below the present 
experimental upper limit, we find that BR{t jj^y) < 10~^, apart from very particular 
parameter points, which may be of special theoretical interest. Only in these cases might 
BR{t fi'j) be within the reach of present experiments at B factories and the LHC. 
The ratios BR{ii — ^ ij'j)/ BR{£k — ^ £ml) depend, within our assumptions, on only one 
unknown complex parameter that appears in the neutrino Yukawa couplings. If possi- 
ble, a measurement of BR{t ^'y)/BR{fi 67) would provide particularly interesting 
information on complex structure in the seesaw model. 



2 Decoupling Assumptions for Yukawa Couplings 

In the seesaw mechanism the observed neutrino masses and mixing are determined by the 
neutrino Yukawa couplings Y^'^ and the mass matrix M/j of the heavy singlet neutrinos. 
The mass matrix of the light left-chiral neutrinos is given by: 



where {H) is the electroweak symmetry breaking vacuum expectation value of the relevant 
Higgs doublet and the complex symmetric matrix C is the coefficient of the dimension-five 
operator CeS = —Cab{LaH){LbH) /2 resulting from integrating out heavy right-chiral 
neutrinos 0. The matrix C, which can be diagonalized by a unitary matrix U^CU^ = 
diag(c^, C2, C3), is at the high scale given by 



(1) 



C 



u 



(2) 
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Similarly as for the quark sector, the number of observables in the light neutrino sector 
is smaller than the number of free parameters in Yj, and Mr. Indeed, even in the basis 
in which M/j = diag(Mi, M2, M3) and Yg = diag(|/e, |/^, ?/t), the most general solution of 
the relation (l2l) for Y„ reads: 



Yt^ = ^MTY.^ADCo^T ■ (3) 

c 

A complex orthogonal matrix accounts for the six-parameter ambiguity in translating 
the 3 neutrino masses and 6 parameters in into the 15 parameters of the neutrino 
Yukawa coupling [g]. Thus, even assuming that the elements of the mixing matrix 
\Ju are determined with sufficient accuracy, the predictions of the CMSSM for flavour- 
violating processes in the lepton sector depend on several unknown factors: the pattern 
of the light neutrino masses (hierarchical, inversely hierarchical or degenerate), the right- 
chiral neutrino masses Ma and on the matrix f2. In practice, one should also remember 
that the Majorana phases in U^, may never be measured. 

In a large number of papers, the predictions for the charged lepton decays have been 
discussed under various specific assumptions about those unknown factors [H 13 El 13 UHl 
EHlEllESEniElEniEElEIlEHlEni. in the present paper we re-examine the predictions 
for the charged lepton decays of the CMSSM under two very general assumptions. We 
assume the hierarchy in mass for both, left- and right-chiral neutrinos 

Ml < M2 < M3 (4) 



We take m^g = yAmatm = 0.05 eV, rrij,^ = ^ Am'^^i = 0.008 eV and the neutrino mixing 
matrix in the form: 



C12 S12 sise 



-7i + --- + 73 I -diagle'^Se^^M) (5) 

\ ^ _L _£12. I J_ 



where = sin^ 6^12 = 0.315, the dots stand for small terms ~ O^sis) (for the sake of 
definiteness we have assumed that < 6*23 < vr/2). For the U;',^ entry we use the two 
scenarios specified in Table 1. We do not assume any particular values of the Majorana 
phases 0i and 02, but treat them as free parameters. It is important to note that, for 
hierarchical light neutrino masses, renormalization effects on the mass eigenvalues and on 
the mixing angles, in particular, are small and can be neglected (see e.g. [^01 El])- Thus 
one can use the measured in Q. 

It has frequently been observed that the simultaneous appearance of hierarchical light 
neutrino masses and two large mixing angles is not 'natural' in the seesaw mechanism 
(see e.g. ^5)- Important exceptions are so-called sequential neutrino dominance models 
with one neutrino decoupled ^Hl or even absent [321 EBl ISEl El]- Otherwise some 
particularities, as found, e.g., in models with horizontal symmetries [33 EEl EI|, are 
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Table 1: Selected values of\J]^ at the electroweak scale used in the analysis. 

necessary. This motivates our second assumption that one singlet right-chiral neutrino, 
not necessarily the heaviest one, decouples from the seesaw mechanism, in the sense that 
at most one neutrino, say Na, contributes to the mass rriu^ of the lightest neutrino. As 
has been explained in [HHlEni, for Y^'^/Ma and Y^/Ma the matrix ^^ in (jSI) 
has to be, depending on the index A, of one of the following three forms: 

/I \ 

decoupling oiNiQ,= \Qzp\ (6) 

V TP ±z J 

/ z p \ 

decoupling of = 1 (7) 



decoupling of A^s ft = \ ±p \ (8) 

Woo; 

where T p'^ = 1. Since ^Iab determines directly the contribution of the right-chiral 
neutrino Na to the mass of the light neutrino ub [2l|i the assumption that Na effectively 
decouples means that it is only Na that contributes to m,^^ , and that it does not contribute 
to m^jg. In the limit of strict decoupling, in which Y-^^/Ma — for i? = 1,2,3, the 
lightest left-chiral neutrino is massless because, as follows from ©-([HI) and Q, ^JrfvJ[ oc 
y\5i)^^ I \/ Ma- Of course, in realistic cases the zeroes in O are non-zero numbers ^ 1. 
The stability of the patterns Q-JHl) with respect to radiative corrections is analyzed in 
the Appendix. 

Another motivation for seesaw models with the light-neutrino masses dominated by 
two heavy singlet neutrinos follows from the hypothesis that the cosmological inflaton 
field is one of the heavy sneutrinos, Na-, say. As has been discussed in jT2|, in order 
to reproduce the measured characteristics of the CMB in such a scenario, and to agree 
with the data on the cosmological structure formation, the inflaton-sneutrino and its 
superpartner, the right-chiral neutrino N a must have a mass Ma — 2 x 10^^ GeV. The 
reheating temperature following inflaton decay is 

■Ma 



Trh ~ JruAMp—^, (9) 



where 
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mA ^ (Y.Yt) {H)yMA = E l^ABl'm,, - (10) 

B=l 
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For Ma — 2x 10^'^ GeV, it follows from ^ that one needs ~ 10~^^ eV in order to obtain 
Trh ~ 10^ GeV 1101. This in turn enforces a pattern in the Q. matrix with (1, 0, 0) in the 
y4'th row. Thus reheating temperature Tp^n after inflation that is low enough to solve the 
cosmological gravitino problem is possible in sneutrino inflation scenarios provided the 
neutrino partner of the Na inflaton (not necessarily the heaviest one) decouples from the 
seesaw mechanism. The possibility of explaining the baryon asymmetry of the Universe 
through leptogenesis following inflation with the patterns ©-(jHI) will be discussed in Sec. 

m 

The three patterns of Vt shown in Q, JH) and ^ provide distinctly different realiza- 
tions of the seesaw mechanism. However, as we explain later, they give similar predictions 
for LFV decays in a large range of parameter space. It follows that the predictions for neu- 
trino masses and mixings and for LFV are invariant under certain transformations on the 
neutrino Yukawa matrix. Therefore, even with all possible data on neutrino masses and 
mixings, as well as on LFV processes that can (in principle) be obtained experimentally, 
the bottom-up approach can determine Y^, only up to these transformations. 

The analysis presented in this paper is based only on the general assumptions described 
above. It is, however, worth commenting on the link of such a bottom-up approach to 
model-dependent (top-down) studies based on speciflc 'theories' of the neutrino Yukawa 
matrix and of the masses of the heavy singlet neutrinos, based for example on texture 
zeroes. Once such a theory of flavour is specified in some electroweak basis, as in [23 IIE] 
for example, Q. can be read off easily from (jSJ, by first performing the necessary field 
rotations diagonalizing the heavy singlet neutrino mass matrix and the charged lepton 
Yukawa coupling matrix. One may envisage a situation in which the electroweak basis in 
which the fiavour theory predicts texture zeroes in Y^^ is only slightly different from the 
basis in which our assumptions are formulated, i.e., the matrices Yg and Mr are almost 
diagonal. The equation Q fixes then the value of the parameter z in terms of the ratios 
of some elements of and of 171^2 /mu^. Indeed, for the pattern (P) [pattern 
Yl^ [Y^^, Yl^] can be neglected, and if in addition the texture gives Yl^ = [Y^-^ = 0] 
for 5 = 2,3 then one obtains from Q 

±- = l-i^(^^l - 1 T 0.06e--^^ . (11) 




Similarly, if Yl^ = [Y^^ = 0, Yl^ = 0] for B = 2,3 because of texture zeroes, then 

±0.3e-''^^ . (12) 
For a diagonal matrix M/j, the solutions (fT^ and (fTTll remain approximately valid if 



lUf^UJ^ for A<B, (13) 

where is the matrix diagonalizing YlY^. The solution (fTT|) is more sensitive to de- 
partures of Nlji from the diagonal form than is (fT2|l . It will be useful to remember these 
particular qualitative patterns with z ^ 0.3 or z ~ 1 in our subsequent analysis. 
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3 Implications of Decoupling for LFV Decays 



In this Section we first present the general background for calculating the branching ratios 
for the LFV decays Ia ^ IbI and then discuss in detail the predictions that follow from 
the pattern Q of the matrix f2. The predictions following from the two other patterns 
are similar and will require only a short discussion. For each of the patterns we also 
discuss possible scenarios for leptogenesis. 



3.1 General Formalism 

The branching ratios for LFV decays are well described by a single-mass-insertion ap- 
proximation |2nl l2i] : 

BR{Ia ^ Ibi) ~ 7^-^("^o, Ml/2, /i) |mi^5 1' tan^ /3, (14) 

where is a function of the soft supersymmetry-breaking masses fixed at the high scale 
Mx. The ratio /(mo,Mi/2) = J^(mo, Mi/2)/J^(100 GeV, 500 GeV) for tan/5 = 10 and 
/i > is shown Fig. [T^. The reference values mo = 100 GeV, M1/2 = 500 GeV are 
chosen for their consistency with the astrophysical cold dark matter constraint [11^. In 
the subsequent figures, we show predictions with the factor /(mo,Mi/2) removed, and 
Fig. can be used to obtain the corresponding predictions for other values of mo,Mi/2. 

The off-diagonal entries in the slepton mass matrix m| are generated radiatively by 
the renormalization-group (RG) evolution from Mx to the electroweak scale. In the 
commonly-used approximation of using only a single recurrence in solving the relevant 
RGBs, they are related to the neutrino Yukawa couplings by 



m? .„ = «:^(Y,^^)*(At + A£c;)Y,^^, (15) 



'■Lab 

C 

where k = -6ml " 2^o, = ln(Mx/M3)/167r2, Aic = IniMs/Mc) /16tt^ i42j. This 
result is not very accurate, as we demonstrate in Fig. [T]b, where we show the ratio of (fTH| 
to rhi^^ as obtained by solving numerically the one-loop RGEs. Although the results 
presented in this paper are based on the exact numerical integration of the RGEs with 
successive decoupling of the singlet neutrinos at the appropriate scales, we note that the 
numerical result is well approximated by the following correction to k: 



36 

- — (?2(2M2/2 + 3ml + ^0 + 2A0M1/2) 
+yK3Qrnl + 2AAl) + {Qml + 12Al)Tr(YlY^)] At, (16) 



where the top-quark Yukawa yt and the other couplings are taken at the scale Mx- The 
correction (fT6|l is obtained by including the second recurrence in solving the RGEs. How- 
ever, the approximation (|THll is sufficient for the qualitative discussion of our results. 
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Combining (fTK|l with ^ and (|H) and using one of the decoupling scenarios (jHI)-®, one 
may obtain formulae for the decay rates BR{Ia Ibi)- It is useful to discuss first the 
results obtained under the technical assumption that the hierarchies of neutrino masses 
obey 

m^^M2 < m^^M^ for pattern (jHI), 

rriy^Mi < m^^Ms for pattern 0, or (17) 

rriy^Mi < my^M2 for pattern (jH]), 

where the inequality signifies a ratio of at least a factor 2 or 3. We study deviations from 
(fT7|l at the end of Subsection As follows from eq. under these assumptions and 
for generic values of the z and p parameters, the formula (fTHll is well approximated by the 
contribution of the product (Yj^^)*Yi^^ for the structures © and 0, and of the product 
{Y'l'^yY'l^ for the structure ((HI). Nevertheless, for the sake of the future discussion, in 
the formulae given below we display also the nonleading contribution to the formula (fT5|l . 
In our numerical calculation we keep, of course, complete mi^^ obtained by solving the 
RGBs. 



3.2 The Decoupling of A^i and LFV Decays 



With the pattern (P), the mass insertion ml^^ relevant for r — > /i7 reads: 

'\jf\]f* [\z\^ + S\l- z'\) +R\Jf\]f* [SzVl^* - z*VT^) + 



mL - K^^^^x (18) 



+RVfVf* {Sz*VT^ - zVl^*) + VfUf* {S\z\' + |1 



-z' 



where R = ^Jm^Jm^^ ~ 0.41 and S = M2(l+A£2/At)/M3 ~ M2/M3 gives the subleading 
contribution of the product (Yl'^)*Yl^ . The branching ratio does not depend strongly on 
the masses Mi and M2 and the Majorana phase 0i, as long as Mi < M2 <C M3, i.e., for 
S* < 1. For the results presented below, we take M3 = 5 x lO^'' GeV, M2 = 3 x 10^^ GeV 
{S ~ 0.1) and Mi = 2 x 10^^ GeV, consistent with infiation being driven by the lightest 
singlet sneutrino. Our choice of M3 makes Yj^^ naturally of order of unity. Moreover, 
due to maximal atmospheric neutrino mixing, Re(U|^^U^^* — Ul^Ul^*) ~ 0(si3) and the 
real parts of the second and third term in (fT8|l cancel each other almost completely. Then 
it is only \m{rnj^^^) which depends on the Majorana phase 02, and the branching ratio 
BR{t — > /i7) oc |m|^2p consist of a sum of two positive terms: a constant one and a term 
which oscillates with 02 as sin^[02 — arg {z*^/l — z"^)]. The relative magnitude of these two 
terms depends on the phase of z, and a non-zero value of U^^ adds a slight modulation 
proportional to — sin (02 — arg {z* ^/l^^z^)) . 

This behaviour is clearly seen in Fig. [2l which shows BR{t ^•y) as a function of 
02 for |zp = 1/2, tan/5 = 10 and Aq = 0, for the two different values of the U^^ element 
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specified in Table [H As already mentioned, the values of the branching ratio for other 
values of mo and M1/2 can be obtained by appropriate rescalings that can be read off from 
Fig.ll^. The dependence of BR{r — > fi^y) on \z\ is depicted in Fig. [HI For each value of \z\, 
the maximal and minimal predictions for this branching ratio are shown. In particular, 
the allowed range for = 1/2 corresponds to a horizontal squeezing of Fig. El 

At this point, it is worth summarizing the dependence of BR{t fij) on the model 
parameters. The dependence on \z\ gives a factor that varies by 10^ for \z\ changing 
between 0.01 and 1.4. As can be seen in Fig. UK, for mo and M1/2 varying in the range 
100 GeV— >1 TeV, the decay rate can change by another a factor of roughly 10^. Further- 
more, the dependence of (fT^ on tan^/3 brings in also a factor of order a few to ~ 10^. 
One should also remember the quadratic dependence of BR{t — » /i7) on M3. The fact 
that 02 is inaccessible to experiment introduces only a small additional uncertainty of the 
order of factor 4. 

Collecting all this information, we see that, under the present assumptions and with 
the anticipated experimental sensitivity down to BR{t — > yU7) ~ 10~^, this decay can in 
the scenario © be experimentally accessible for a large range of parameters. The question 
which part of this parameter space can be consistent with the present experimental bound 
BR{fi 67) < 1.2 X 10~^^ is considered below. It is however worth noting that, assuming 
that the soft masses and tan /3 can be determined by other measurements with sufficient 
precision, the measurement of this decay would, owing to the strong dependence on \z\: 
BR{t fi'-y) oc \z\*, provide a relatively narrow range of \z\ for fixed masses of the heavy 
singlet neutrinos. 

The quantities m^j^ relevant for /i — > 67 and mj^.^-^ relevant for r — »• 67 can be 
approximated by: 



_ 2 m^,^M3At 



RUfUl'* [SzVT^* - z*VT^) 



'Lai - - 

+R^\jfvl'* {S\z\' +\l-z^\)+ Vfvl^* {\z\' + S\l- z^\)] , (19) 

where we have dropped the terms suppressed by both U^^ and R. For the sake of qual- 
itative discussion, the third term can always be neglected (for U^^ 7^ 0, this is a good 
approximation provided |1 — is not too close to zero). For BR{fi — > 67), the factors de- 
pending on the neutrino masses and mixing in the first two terms in (fT9| can be estimated 
as: 

i^lU^^U^^I ~ 0.16 i?^|UfU^2| ~ 0.05, (20) 

whose ratio is 3.1 ± 0.5. For S* <^ 1 the relative phase between the first two terms in 
(fT9ll is 02 — arg (-2*Vl — 2^) + tt. This is clearly seen in Figs. andlUb, where we show 
BR(n — i> 67) //(mo, M1/2) as a function of the Majorana phase 02 for = 1/a/2 and three 
representative values of arg^;. The strength of the 02-dependent interference pattern seen 
in Fig. m varies with the values of \z\ and arg^;. In particular, it decreases both for \z\ — > 
and 1 2; I ^ 1, where the second term in (fT9|l dominates over the first one. Comparison of 
Fig. Eti with 131b shows that a non-zero value of U;',^ enhances the interference pattern but 
does not change it qualitatively. 
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The dependence of BR{^ 67) on the other parameters is similar to that of BR{t — > 
/i7). In Fig. Owe show BR{fi — > 67) as a function of \z\ for the same values of the other 
parameters as we used in Fig. 121 for BR{t fi'y). The comparison of Figs. [T] and [HI 
shows that, except for two special points at \z\ ~ 0.3 and \z\ ~ 1, the branching ratio of 
the decay /i — 67 for /(mo,Mi/2) > 1 and tan/? > 10 is above the experimental upper 
bound 1.2 x 10~^^. The minimum of BR{jj — > 67) at \z\ ~ 0.3 is due to the destructive 
interference of the first two terms in (fT9|l : for S* -C 1 they are of equal magnitude just 
for 1^1 ~ 0.3, and for any value of arg2; the Majorana phase 02 can be chosen so that 
these two terms approximately cancel each other. The minima for \z\ ~ 1 are instead 
caused by the simultaneous vanishing of the two first terms in the limit S = 0. Therefore 
they appear only for arg2; ^ and vr. For U^^ = the non-zero value of BR{jjL 67) 
at the single minimum seen in Fig. Eti is due to S* 7^ in the second term in (fT9|) : for 
7^ (Fig- ISb) there are two minima whose position is determined by the cancellation 
of the first and third terms in (fT9|) . which holds for z ^ ±(1 — iQ;exp(2i(5 — 02)), where 
a = \\]f*/R\]f*\ ^ 1/5, i.e. for \z\ ^ 1 - |a cos(2(5 - 02))- The minima of Si? (/i ^ 67) 
occur at the values of z corresponding to the texture zeroes (fTT|l and (fT^ . This is obvious 
for Y^^ = 0, since BR{Ia ^ hi) oc jY^^^Y^^p, but it is a non-trivial result for Y^^ = 0, 
which follows from the structure of the masses and mixings of the light neutrinos. 

As follows from Fig.HJ for BR{fi ej) to be consistent with the experimental bound 
in a larger range of \z\ values, the gaugino mass M1/2 at the high scale should be bigger 
than 500 GeV and/or M3 significantly smaller than 5 x 10^"^ GeV. On the other hand, 
Ml/2 ~ 500 GeV corresponds to third-generation masses above 1 TeV at the electroweak 
scale, which begins to confiict with the naturalness requirement. This confiict becomes 
even sharper for higher values of tan/3. This discussion suggests that, if fl is of the form 
the rate of — > 67 decay should be close to the present experimental bound. 

At this point, we can also address the question whether the present experimental bound 
BR{^ — > 67) < 1.2 X 10~^^ allows for BR{t /X7) > 10~^. This can best be discussed 
by studying the predictions for the ratio BR{n —>■ e'j)/BR{T yU7), which depends only 
on z and 02, as the dependence on the soft supersymmetry-breaking parameters, tan/? 
and M3 cancel out. Furthermore, since BR{t /i7) is a monotonic function of \z\ and 
can be changed by a factor of 10 at most by changing the Majorana phase 02, the minima 
of the ratio BR{fi e''y)/BR{T — > ^^7) follow the minima of BR{fi 67). Maximal 
and minimal possible values of BR{jj, e-^)/ BR{t ^■y) obtained by varying 02 in the 
range (0, 27r) are shown in Fig. IHlas a function of \z\ for three different values of argz. 

As seen in Fig. IHl generically BR{ji e'y)/BR{T — > ^^7) ~ 0.1 to 1, that is BR{t — > 
/i7) ^ 10~^° for BR{fi — ^ 67) ^ 10^^^. However, there are exceptions to this rule in a 
few isolated regions corresponding to the minima of BR{ji — > 67) seen in Fig. This 
means that observation of r — /i7 at a rate > 10^ in future experiments would be a strong 
constraint for top-down models of neutrino masses and mixings. 

Finally, we briefiy consider the decay r — > 67. The dependences of the predictions for 
BR{t — > 67) on the soft mass parameters z and 02 are similar to those of BR{n 67). 
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However, since the relative phase between U^^ and Ul^ approximately equals vr, for the 
value of 02 for which BR{jj, — > 67) is minimized BR{t 67) is maximal and vice versa, 
as shown in Fig. The dependence of BR{t e'y)/BR{T — > fi'j) is depicted in Fig. |H1 

We have presented predictions obtained under the supplementary assumption (|T7|l . 
However, for M2 close enough to M3 but still reasonably smaller, the first condition (fTTjl is 
no longer satisfied. It is therefore worthwhile to check the behaviour of the predictions for 
the decays rates when M2 approaches M3. In particular, one may wonder whether the deep 
minima in BR{jjL 67) at \z\ ~ 0.3 and \z\ ~ 1 are filled in by additional contributions 
to (fTKll coming from the product (Y^^)*Y^^, which may now be non-negligible. This is 
shown in Fig.lHl The two minima persist even for relatively large values of the mass ratio, 
M2/M3 < 0.4 for = and M2/M3 < 0.25 for \\Jf\ = 0.1, but the first minimum shifts 
to slightly higher values of \z\ for larger values of M2/M3. 

In this Section we have discussed the results for BR{Ia — ^ Ibi) obtained under the 
hypothesis that the lightest right-chiral neutrino decouples from the seesaw mechanism. 
To a very good approximation, these results do not depend on the mass Mi. The particular 
value Ml = 2 x 10^^ GeV, is compatible with infiation driven by the lightest sneutrino. As 
analyzed in jTZ], succesful nonthermal leptogenesis with low reheating temperature can 
then occur. This is because for the zeroes in © representing small complex numbers it is 
easy to obtain values of the CP asymmetry parameter ei saturating the Davidson-Ibarra 
upper bound [i3|^. 

For any other value. Mi 7^ 2 x 10^^ GeV but still with the pattern some other field 
must be responsible for infiation. Disregarding the cosmological gravitino problem, one 
may then contemplate the possibility of conventional thermal leptogenesis. It turns out 
that, with zeroes in Q representing small complex numbers, leptogenesis can be realized 
for Trh ~ 10^ GeV In both cases, the final lepton number asymmetry does not 

depend on the parameter z of the fi matrix. The possibility of lowering the reheating 
temperature necessary for the thermal leptogenesis in a scenario with degenerate neutrinos 
A^i (decoupled) and N2 (not decoupled) remains an open question for the pattern 

3.3 LFV decays and the decoupling of A^2 or A^3. 

The predictions for the LFV decays with the decoupling of A^'2 [pattern 0] or A^3 [pattern 
((HI)] are easy to discuss if one remembers that mi^^ in eq. (fT5|l is dominated by the 
contribution of the product {Y^'^yY^^ or (Y^^)*Y^^, respectively. This means that in 
the case ^ the predictions for the decay rates are the same as for pattern for the same 
values of Yj^^ and if in both cases the subleading contributions to are negligible ^. 

For the decoupling of iVs, the predictions are the same as for the decoupling of A''i if we 

^However, for moderately hierarchical heavy neutrino masses it is possible to obtain larger values of 
the CP asymmetry than the abovementioned bound jlij . 

^Predictions in the two patterns are the same including the subleading effects if the Y^"^ in the case 
10 are numerically equal to the Y^'^ in the case 
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change Yj^^ — Y^"^ and keep the numerical values unchanged. At this point, it is worth 
recalling equation Q: 



r2,A 

V 

r2A 

V 



iMl''' (l^siClUf * + 032C2Uf * + ^^33C3U 
iMl'^ (f^2lClUf * + 022C2Uf * + f^23C3U 



A3* 

V 

A3* 

V 



(21) 



We see that the predictions following from decoupling of A''i and N2 are, to a good ap- 
proximation, identical for the same values of M3. In particular, the mass of the decoupled 
neutrino can be kept at 2 x 10^^ GeV (for sneutrino inflation) Decoupling of iVs leads 
to the same LFV decay rates as decoupling of A''i for M2 replaced by a numerically equal 



M3. Moreover, with M3 = 2 x 10 GeV the dominant contribution to Y;; in ()2T]l comes 
from M2 < M3. So the BR{Ia Ibi) discussed in Section El are rescaled by the ratio 
Ml /Ml Since, for P, M3 > 2 x GeV and, for ©, M2 < M3 = 2 x 10^^ Q^y^ ^.^e 
predictions for the decay rates are strongly suppressed in the latter case. We conclude 
that, if the heaviest right-chiral neutrino is the inflaton, the decays Ia — IbI are unlikely 
to be observed experimentally. 

Let us now consider the possible types of leptogenesis with N2 or A^3 decoupled. In 
nonthermal leptogenesis following sneutrino-driven inflation |17| . the final lepton asymme- 
try is generated directly in the decay of the vacuum condensate of the sneutrino-inflaton 
field N2 or A''3. Since the flnal lepton- number-to-entropy ratio Yl generated in this way 
is given by: 

Yl = ^eA^ , (22) 



a low reheating temperature Trh 



4 ^Minfl 
10^ GeV and Mnfl 



2 X 10^3 GeV require the CP 
asymmetry parameters €2 or 63, respectively, to assume values of the order of eref/10, 
where eref = ZTh)^' ~ 4 X 10-3 (for iVi-driven inflation eref is the same as the Davidson- 
Ibarra bound on ei |1S1)- In Table El we display the one-loop contributions to eA/^reiVAE 
from the individual virtual states Ne (and Ne) [IE], where 



Vae 



ED) 



<E 



m 



I'D 



n 



D ^1^3 



ED\ 



(23) 



For the pattern only the contribution from Ni is suppressed by the small value of 
{M1/M2Y, whereas for the pattern ^ all the contributions are small. This suppression 
can only be overcome if tjae, that is some iri^;/)^, are large, which, recalling the interpre- 
tation of Q as the dominance matrix, means some flne-tuning in the seesaw mechanism. 
This can be seen in Figure UHl where we show in the form of a scatter plot the CP asym- 
metries 62 and €3 generated in the direct decays of the sneutrino-inflaton N2 [Ns], relevant 
for the patterns ^ respectively, as functions of the parameter p = \/\ — of the 

matrix. We have chosen argp = 7r/4 and argf2B2 = ai'gf^_B3 = 7r/4 for B = 2,3 and 
scanned over |r2^2| and \^b3\ in the range 10-^ — 10^^ with a uniform distribution over the 



^With the exception that, in the case (O, subleading effects may become important only at the expense 
of lowering M3, i.e., when the LFV decay rates are suppressed. 
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pattern ^ 


pattern ^ 


E = 1 
E = 2 
E = 3 


1 


3M| V Afs J 

H4fl + ln^l 





Table 2: Estimated maximal contributions to \e2 / ereiil2E\ for the pattern and to 
\^z/^vciV^E\ for the pattern due to virtual Ne exchange in the loop (see the text for 
explanation) . 

logarithmic scale. It is clear from Figures ITTIk and ITOb that, whilst the Davidson-Ibarra 
bound can be saturated for the pattern with l-z], |p| ^ 1, for the pattern ((HI) this 
is possible only for |p| > 3 and, hence, requires fine tuning in the seesaw mechanism. 

Finally, we can abandon the hypothesis of sneutrino-driven infiation and consider more 
conventional thermal leptogenesis with the patterns and (jH)). The masses of all the 
three right-chiral neutrinos, including the decoupled one, can be then arbitrary. If they 
are split sufficiently for the final lepton number asymmetry to be generated entirely in the 
decays of the lightest right-chiral neutrino, the only relevant CP asymmetry parameter is 
ei, which in both cases and (jHl), takes the form 



3Mi 


Im(z2) (ml -ml) 




87r(i7)2 


z 









(24) 



The important wash-out parameter mi then reads 

rfii ^ {\z\'^m„^ + \1 — z'^\m^^) . (25) 

The basic consequences of and (f^ were explored in |3H]-^ It was found that, because 
fhi is bounded from below by m^^j and not by m^-^ , the wash-out processes are very efficient 
and successful leptogenesis requires compensation by a large ei, i.e., large Mi. For the 
patterns (0) and ©, with Mi = 10^^ GeV and Mi = 5 x 10" GeV, regions of the plane 
(|2;|,argz) leading to successful thermal leptogenesis are shown in Fig. [TT] ^. The two 
panels of Fig. fTD correspond to the two cases specified in Table 1. Contours of constant 
BR{fj, 67) minimized with respect to 02 are also shown. As can be seen, regions of the 
(|2;|,arg2;) plane where the thermal leptogenesis reproduces the observed baryon number 
asymmetry of the Universe and regions where the experimental bound on BR{jjL — > 67) 
can be respected overlap with each other. However, since Trh > Mi > 10^^ GeV for this 
mechanism to work, the cosmological gravitino problem has to solved in some other way, 
perhaps by the gravitino being the lightest supersymmetric particle (see, however, 

m- 

"In only the decoupling of the heaviest right-chiral neutrino, i.e., the pattern (jHJ, was considered 
explicitly. However, it is clear that these results apply to the pattern 10 as well. 

^Recent advances in the calculation of the lepton number asymmetry ^3 have been included. 
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As has been discussed in |3H1 HHl HHl EDI EH E21 , for the pattern © Trh < 10^ GeV 

can be obtained if Ni and are tightly degenerate, with 1 — M1/M2 ~ 10~^. In this case, 
however, since is decoupled and M2 < 10® GeV, branching ratios BR{li Ij'j) are 
~ 10"^*^, too small to be accessible to current experiments. For the pattern one of the 
two degenerate neutrinos would be the decoupled one, and to obtain succesful leptogenesis 
with TptH ~ 10® GeV one would probably need the degeneracy of all the three right-chiral 
neutrinos. In this case the LFV decays would again be inaccessible. 



4 Conclusions 

We have calculated the branching ratios for the decays Ia ^ IbI in the MSSM under the 
assumption that the masses of both the light and heavy singlet neutrinos are hierarchical, 
and that one singlet neutrino (not necessarily the heaviest one) decouples from the seesaw 
mechanism. For each of the three possible realizations of the decoupling hypothesis we 
have also discussed possible leptogenesis scenarios. 

The predictions for BR{Ia —>■ Ibi) do not depend on which neutrino decouples. To 
a very good approximation they depend on only one non-decoupled neutrino mass, the 
heavier one. Apart from the dependence on the soft mass terms and on tan/3, parame- 
ters that will hopefully be measured in other experiments, the branching ratios depend 
strongly on just one complex parameter 2;, which fixes all the relevant Yukawa couplings. 
Consequently, our two assumptions lead to relatively rigid predictions for the decay rates. 
For mo, M1/2 ~ 1 TeV and M4 > 10^^ GeV, the decay rate for /i ^ 67 is predicted to 
be close to the present experimental bound 1.2 x 10~^^, namely BR{fi — > 67) ^ 10"^''. 
However, for BR{^ e-y) < 1.2 x 10~^^, the decay rate for r /i7 is generically below 
the anticipated experimental sensitivity 10~^, except for some special values of 1^1 . The 
ratios BR{fi — > e'~^)/ BR{t — >• /i7) and BR{t — >• e'y)/BR{T — >• fi'j) depend on z only, with 
the dependences on the other parameters cancelling out. These ratios would therefore be 
particularly interesting to compare with experiment. 
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Appendix - Stability of the fl patterns 



In bottom-up evolution, the Yukawa couplings Y^^, Y^^ and Yl^ enter into the RGEs 
for the soft slepton masses at the scales Mi, M2 and M3, respectively. We see from © 
that they are expressed at these respective scales in terms of the matrix f2, which is 
taken to have one of the scale-independent forms ©-(jHI)- However, once the Y^^ (say) is 
fixed by ^ in terms of ft at the scale Mi, their scale dependences are controlled by the 
appropriate RGEs, and at the scale M2 they do not satisfy the relation ^ with scale- 
independent fl. Thus, fixing Y^^ at M2 from ^ with a scale-independent orthogonal 
matrix fl is not a fully consistent procedure, as it distorts the relations between Yl^ and 
Yl^ which follow from the assumption that, e.g., with ^ Yl^ /Mi ^ for -B = 2, 3. 



We can get some control over the quality of this approximation in the following way. 
Suppose that all three right-chiral neutrinos are integrated out at a scale Qq chosen in 
such a way that the threshold corrections to the relation Q vanish. It is easy to verify 
that in a one-loop approximation Mi < Qq < M3. We denote by Co the value at Qo of the 
coefficient of the dimension-five operator C5 = —Cab{LaH){LbH)/2 related to C(M^) 
directly by the RGEs of the MSSM derived in At the scale Qq, we can use © to fix 
the Yukawa couplings in terms of ft{z) given by one of the three patterns ©-JHI)- Once 
^u{Qo) is determined by fl, the RGEs for Y^, Ma and m^^ can be used to calculate 
these quantities at any other scale Mi < Q < Mgut- Indeed, one can obviously define 
the running of the coefficient of the O5 operator also for Q ^ Qo, because this operator 
could be added to the initial Lagrangian with right-chiral neutrinos present. Above Qo, 
the RGE for C is 



4c = -KC 

dt 



YlY^f + (Yt Y.)^l C - C f(YtYe) + (Yt Y.)! , (26) 



where K = —Qgl — 2gy + 2Tr(3Yj^Y„ + Y^Y^,). By evolving C to an arbitrary scale Q 
by using the RGE with Co as the initial condition at Qq ^ one can use (jH)) to define 
flAB at any scale Q 7^ Qq hy the formula 

Qab = -^(O^Y.U,)^^-, (27) 

y/ Ma Cb 

where \5i, diagonalizes C at Q, and Or is the unitary matrix diagonalizing Mr at Q. The 
question of consistency is now how much (Iab differs from Q-ab- In particular, VtAB{Q) 
defined in this way needs not be orthogonal, because in general Q{Q) would be numerically 
different from YJM^^Yj, at the same scale. 

However, it is easy to see that in supersymmetry 0,ab{Q) defined in this way is 
orthogonal because the RGE (f^ is exactly the same as the RGE of the combination 
YjM^^Yy, as obtained by using the chain differentiation rule and the RGEs of Yj, 

^For hierarchical hght neutrino masses, the renormaUzation effects between the scale Mz and Qo on 
the mass eigenvalues and, in particular, on the mixing matrix are small, and Co can be identified 
with U*diag(TOi^j , TO^2 , 771^3 )Uj, at the electroweak scale. 
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and given, e.g., in pTj. This is because, in supersymmetry, the Yukawa couplings, 
the right-chiral neutrino mass term and the dimension-five operator are all F terms. 
Therefore, due to the non-renormalization theorems, the RG running of their coefficients 
are entirely given by the wave function renormalization of the (super)fields out of which 
they are built. It then follows that renormalization of the Na (super)fields, which does 
not enter the running of C, also cancels out in the expression for the RGE for YjM^^Y,^. 

The RGE for Q can be then obtained by simply differentiating the formula (f27jl . using 
the known RGEs for Ui, and cb [SOIISI], the RGE for (see, e.g., the Appendix of [SI] ^) 
and the RGEs for Ma and Or obtained by applying the technique of [SH EOl IHI] to the 
RGE for Mr derived in One then gets 



^Ma = 4Ma [oiiY^YlfOn 



dt 



AA 



c 



AC^CB 
R 



(28) 



where Eq = —£oi ^o^ = and 



.AB 
'O 



^Ma + Mb^ 
Ma -Mb 



AB Ma 
-2i- ^ 



Mf 



Ma + Mb 



-Im 



for A B. Combining all the elements, one gets 



OUY.YlfOn 



AB 



(29) 



-r^^AB 

dt 



Ma 



EE 

Dj^A C 



Mom, 



+2 



{Mn - Ma) 
Md 



Re (n, 



AC^*DC 



+ i 



Mom^ 



Cj^B D 



+2 ^ Y.^Acy, 

C^B D 



AC 



2 



m 



(m 



m 



-Re 



Qdb + i 



{Md + Ma) 
M, 



n 



m 



Im (ytAC^*DC 

Im {p.*jjc^DB 



DB 



Re(u^ 



DC*-^DB 



+ i 



— Im 



DC*jjDB 



(30) 



where the y1 are the eigenvalues of Y^Yg 



For non-degenerate heavy neutrinos, it is easy to see that for indices AB corresponding 
to the small entries of in a given pattern (e.g., VtiA and Q,ai with A = 2, 3 for the pattern 
(P)), the first two lines are small compared to 1 because they are always suppressed by the 
mass of the lightest left-chiral neutrino rriy^ . The only potentially dangerous contribution 
comes from Sl^c = 1 in the last line. But this would require that C = 1, and this term is 
suppressed at least by a factor ~ ^Jm,J^/m,J^y1, which is also small compared to 1. Hence, 

for the patterns ©-(jHI) ~ independently of the choice of Qq, if the heavy neutrinos 
are not degenerate. 

Stability of the patterns ©-JHI) may be more problematic for arbitrarily tight degen- 
eracy of the right-chiral neutrinos. One may, however, ask whether the pattern (jH)) is 

'^In the expression (B.12) of this reference the factor of 3 should multiply YJY^ instead of Y^Y^. 
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destabilized by the minimal degeneracy of Mi and M2 allowing for succesful leptogenesis 
with Ml < 10^ GeV [MilMl- To make this estimate, we recall that with degenerate Mi and 
M2 the final lepton number asymmetry is roughly the same as would be obtained from 
decays of only one heavy neutrino having the effective mass Mi^ = Mi/ (M2/M1 — 1). For 
the pattern decays of a single neutrino produce the right lepton number asymetry if 
its mass is ^ 10^^ GeV [3^. Therefore, with the degeneracy we also need Mf^ > 10^^ GeV. 
Consider now the contribution of the first line of (|3n|l to the derivative of the small el- 
ement fifii, B = 1,2. The most dangerous contribution for A = 1 comes from C = 3 
and D = 2. It is of order {Mi/ {Hf){M2m^J{M2 - M{) ~ {Mfm^J{H)^) x < 1. 
Similarly the (nonresonant) contribution of -D = 3 is suppressed by small fiss. Therefore 
the pattern (jHl) is not destabilized in this case, either. 

References 

[1] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81 (1998) 
1562; Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 89 (2002) 011301; 
K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90 (2003) 021802 

[2] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57 (1986) 961; L.J. Hall, V.A. Kost- 
elecky and S. Raby, Nucl. Phys. B267 (1986) 415. 

[3] M. Gell-Mann, P. Ramond and R. Slansky, in Proceedings of the Supergravity Stony 
Brook Workshop, New York, 1979 (eds. P. van Nieuvenhuizen and D.Z. Freedman, 
North-Holland, Amsterdam); T. Yanagida, in Proceedings of the Workshop on Uni- 
fied Theories and Baryon Number in the Universe, Tsukuba, Japan, 1979 (eds. A. 
Sawada and A. Sugamoto, KEK Report No. 79-18, Tsukuba); R. Mohapatra and 
G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[4] M. Fukugita and T. Yanagida, Phys. Lett. B174 (1986) 45; M. A. Luty, Phys. Rev. 
D45 (1992) 455. 

[5] W. Buchmiiller, P. Di Bari and M. Pliimacher, Nucl. Phys. B643 (2002) 367, Phys. 
Lett. B547 (2002) 128, Nucl. Phys. B665 (2003) 445 and preprint DESY-030100, 
UAB-FT-551, CERN-TH-2003-199 [arXiv:hep-ph/0401240j 

[6] J. A. Casas and A. Ibarra, Nucl.Phys. B618 (2001) 171 

[7] S. Davidson and A. Ibarra, JHEP 0109 (2001) 013; J. R. Ellis, J. Hisano, M. Raidal 
and Y. Shimizu, Phys. Rev. D66 (2002) 115013 

[8] G.C. Branco, T. Morozumi, B.M. Nobre, M.N. Rebelo, Nucl. Phys. B617 (2001) 475; 
A.S. Joshipura, E. A. Paschos, W. Rodejohann, JHEP 0108 (2001) 029; D. Falcone, 
Phys. Rev. D66 (2002) 053001; G. C. Branco, R. Gonzalez Felipe, F.R. Joaquim, 
M.N. Rebelo, Nucl Phys. B640 (2002) 202; W. Rodejohann, Phys. Lett. B542 
(2002) 100; Eur. Phys. J. C32 (2004) 235; E. Kh. Akhmedov, M. Frigerio, A. Yu. 



16 



Smirnov, JHEP 0309 (2003) 021; S. Pascoli, S. T. Petcov and W. Rodejohann, Phys. 
Rev. D68 (2003) 093007; S. Davidson and R. Kitano, JHEP 0403 (2004) 020 



[9] J. R. Ellis and M. Raidal, Nucl. Phys. B643 (2002) 229 
[10] S. Davidson, .JHEP 0303 (2003) 037 
[11] S. Davidson and A. Ibarra, Nucl. Phys. B648 (2003) 345 
[12] G. C. Branco et al., Phys. Rev. D67 (2003) 073035. 

[13] A. Y. Smirnov, Phys. Rev. D48 (1993) 3264; S.F. King, Phys. Lett. B439 (1998) 
350 and Nucl. Phys. B562 (1999) 57 and Nucl. Phys. B576 (2000) 85 and .JHEP 
0209 (2002) Oil and Rept. Prog. Phys. 67 (2004) 107; S. Davidson and S.F. King, 
Phys. Lett. B445 (1998) 191; Q. Shafi and Z. Tavarkiladze, Phys. Lett. B451 (1999) 
129; G. Altarelli, F. Feruglio and I. Masina, Phys. Lett. B 472 (2000) 382; S.F. King 
and G.G. Ross, Phys. Lett. B520 (2001) 243. 

[14] For a review see, e.g., G. Altarelli and F. Feruglio, Phys. Rept. 320 (1999) 295 and in 
Neutrino Mass, Springer Tracts in Modern Physics, eds. G. Altarelli and K. Winter 

[15] H. Murayama, H. Suzuki, T. Yanagida and J. Yokoyama, Phys. Rev. Lett. 70 (1993) 
1912; H. Murayama, H. Suzuki, T. Yanagida and J. Yokoyama, Phys. Rev. D50 

(1994) 2356; K. Hamaguchi, H. Murayama and T. Yanagida, Phys. Rev. D65 (2002) 
043512 

[16] J. R. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. B145 (1984) 181; J. R. El- 
lis, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B259 (1985) 175; J. R. Ellis, 
D. V. Nanopoulos, K. A. Olive and S. J. Rey, Astropart. Phys. 4 (1996) 371; 
M. Kawasaki and T. Moroi, Prog. Theor. Phys. 93 (1995) 879; T. Moroi, Ph.D. thesis, 
|arXiv:hep-ph/9503210; M. Bolz, A. Brandenburg and W. Buchmiiller, Nucl. Phys. 
B606 (2001) 518; R. Cyburt, J. R. Ellis, B. D. Fields and K. A. Olive, Phys.Rev. 
D67 (2003) 103521 

[17] J. R. Ellis, M. Raidal and T. Yanagida, Phys. Lett. B581 (2004) 9 

[18] C. L. Bennett et al., Astrophys. J. Suppl. 148 (2003) 1; D. N. Spergel et al, Astro- 
phys. J. Suppl. 148 (2003) 175; H. V. Peiris et al, Astrophys. J. Suppl. 148 (2003) 
213 

[19] G. Lazarides, Q. Shafi, Phys. Lett. B258 (1991) 305. 

[20] J. Hisano, T. Moroi, K. Tobe, M. Yamaguchi and T. Yanagida, Phys. Lett. B357 

(1995) 579; J. Hisano, T. Moroi, K. Tobe and M. Yamaguchi, Phys. Rev. D53 (1996) 
2442; J. Hisano and D. Nomura, Phys. Rev. D59 (1999) 116005. 

[21] X. J. Bi, B. Feng and X. M. Zhang, |arXiv:hep-ph /0309195J . 

[22] T. Blazek and S. F. King, Nucl. Phys. B662 (2003) 359 



17 



[23] J. R. Ellis, J. Hisano, S. Lola and M. Raidal, Nucl. Phys. B621 (2002) 208 

[24] S. Lavignac, I. Masina and C. A. Savoy, Phys. Lett. B520 (2001) 269 and Nucl. 
Phys. B633 (2002) 139 

[25] S. Pascoli, S. T. Petcov and C. E. Yaguna, Phys. Lett. B564 (2003) 241; S. T. Petcov, 
S. Profumo, Y. Takanishi and C. E. Yaguna, Nucl. Phys. B676 (2004) 453 

[26] A. Ibarra and G. G. Ross, preprint CERN-TH-2003-294, |arXiv:hep-ph/0312138j . 

[27] K. Tobe, J. D. Wells and T. Yanagida, Phys. Rev. D69 (2004) 035010 

[28] A. Masiero, S. Profumo, S. K. Vempati and C. E. Yaguna, JHEP 0403 (2004) 046 

[29] B. A. Campbell, D. W. Maybury and B. Murakami, JHEP 0403 (2004) 052 

[30] P. H. Chankowski, W. Krolikowski and S. Pokorski, Phys. Lett. B473 (2000) 109; 
J. A. Casas, J. R. Espinosa, A. Ibarra and I. Navarro, Nucl. Phys. B573 (2000) 652. 

[31] P. H. Chankowski and S. Pokorski, Int. J. Mod. Phys. A17 (2002) 575 

[32] P. H. Frampton, S. L. Glashow and T. Yanagida, Phys. Lett. B548 (2002) 119; 
M. Raidal and A. Strumia, Phys. Lett. B553 (2003) 72; S. F. King, Phys. Rev. D67 
(2003) 113010; S. Raby, Phys. Lett. B561 (2003) 119; B. Dutta and R. N. Mohapatra, 
Phys. Rev. D68 (2003) 056006; V. Barger, D. A. Dicus, H. J. He and T. j. Li, Phys. 
Lett. B583 (2004) 173; W. 1. Guo and Z. z. Xing, Phys. Lett. B583 (2004) 163 

[33] R. Barbieri, T. Hambye and A. Romanino, JHEP 0303 (2003) 017 

[34] R. Kuchimanchi and R. N. Mohapatra, Phys. Rev.B66 (2002) 051301; R. Kuchi- 
manchi and R. N. Mohapatra, Phys. Lett. B552 (2003) 198 

[35] G. Altarelli and F. Feruglio, Phys. Lett. B451 (1999) 388 

[36] Y. Grossman, Y. Nir and Y. Shadmi, JHEP 9810 (1998) 007; Y. Nir and Y. Shadmi, 
JHEP 9905 (1999) 023; J. Feng, Y. Nir and Y. Shadmi, Phys. Rev. D61 (2000) 
113005. 

[37] M. Hirsch and S. F. King, Phys. Rev. D64 (2001) 113005 

[38] P. H. Chankowski and K. Turzyhski, Phys. Lett. B570 (2003) 198 

[39] A. Ibarra and G. G. Ross, Phys. Lett. B575 (2003) 279 

[40] M. Y. Khlopov and A. D. Linde, Phys. Lett. B138 (1984) 265; I. V. Falomkin et al. 
Nuovo Cim. A79 (1984) 193 [Yad. Fiz. 39 (1984) 990]; M. Y. Khlopov, Y. L. Levitan, 
E. V. Sedelnikov and I. M. Sobol, Yad. Fiz. 57 (1994) 1466; K. Kohri Phys. Rev. 
D64 (2 001) 04 3515; M. Kawasaki, K. Kohri and T. Moroi, preprint RESCEU-04-04 
[arXivTastro-ph /0402490| . 



18 



[41] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Lett. B565 (2003) 176; 
preprint CERN-TH-2003-310 | arXiv:hep-ph/0 312262|. 

J. R. Ellis, J. Hisano, M. Raidal and Y. Shimizu, Phys. Lett. B528 (2002) 86 
S. Davidson and A. Ibarra, Phys. Lett. B535 (2002) 25 

T. Hambye, Y. Lin, A. Notari, M. Papucci and A. Strumia, preprint FUP-TH-2003-48 
arXiv:hep-ph/0312203 | 

G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, preprint CERN-TH- 
2003-240 ^arXiv:hep-ph/0310123j . 

L. Covi, E. Roulet and F. Vissani, Phys. Lett. B384 (1996) 169; A. Pilaftsis, Int. J. 
Mod. Phys. A14 (1999) 1811 

M. Fujii, M. Ibe and T. Yanagida, Phys. Lett. B579 (2004) 6 

M. Flanz, E. A. Paschos and U. Sarkar, Phys. Lett. B345 (1995) 248, Erratum - 
ibid. B382 (1996) 447; M. Flanz, E. A. Paschos, U. Sarkar and J. Weiss, Phys. Lett. 
B389 (1996) 693 

A. Pilaftsis and T. E. Underwood, preprint MC-TH-2003-09 |arXiv:liep-ph/0309342j . 



R. Gonzalez Felipe, F. R. Joaquim and B. M. Nobre, arXiv:hep-ph/0311029 



J. R. Ellis, M. Raidal and T. Yanagida, Phys. Lett. B546 (2002) 228. 



K. Turzynski, preprint IFT-04/05 |arXiv:hep-ph/0401219| . 

P. H. Chankowski and Z. Pluciennik, Phys. Lett. B316 (1993) 312; K. S. Babu, 
C. N. Leung and J. Pantaleone, Phys. Lett. B319 (1993) 191. 

K. S. Babu, Z. Phys. C35 (1987) 69. 

J. A. Casas, J. R. Espinosa, A. Ibarra and I. Navarro, Nucl.Phys. B569 (2000) 82. 



19 



200 



400 



800 



0.7 



200 



400 



600 



800 

M,/j(GeV) 



Figure 1: Left panel: /(mo,Mi/2) = J?-'(mo, Mi/2)/^(100 GeV, 500 GeV) as a function 
of Ml/2 for tan/3 = 10 and /i > 0. The normalization is such that /(mo, M1/2) = 1 
for mo = 100 GeV, M1/2 = 500 GeV, consistent with the dark matter constraints UTj . 
Right panel: The ratio of ih\^^{fu\\) obtained by numerical integration of the full set of the 
RGEs to m|^2(approx.) obtained with the use of the formula ( [J .5) J as a function of M1/2 
for tan/3 = 10 and Aq = 0. 
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Figure 2: BR{t /i7)//(mo, M1/2) as a function 0/ 02 for the choice of the remaining 
phases described in Tabled \z\ = tan/3 = 10 and Aq = 0. Dotted, dashed and solid 

lines correspond to a.Tgz = 0,7r/4, 7r/2, respectively. 
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Figure 3: Extremal values of BR{t — > fi'y) / f{mo, M1/2) as a function of\z\ for the choice 
of the parameters described in Taft/eQ tan/3 = 10 and Aq = 0. Dotted, dashed and solid 
lines correspond to aigz = 0,7r/4, 7r/2, respectively. 




Figure 4: BR{fi — > e7)//(mo, M1/2) as a function 0/ 02 for the choice of the remaining 
phases described in Table\^ \z\ = tan/? = 10 and Aq = 0. Dotted, dashed and solid 

lines correspond to aigz = 0, rc/A and vr/2, respectively. 
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Figure 5: Extremal values of BR{jjL — > cy)/ f{mQ, M1/2) as a function of \z\ for the choice 
of the parameters described in Tabled tan (3 = 10 and Aq = 0. Dotted, dashed and solid 
lines correspond to aigz = 0, ir/A and 7r/2, respectively. 




Figure 6: Extremal values of BR{fi — ^ e'y)/BR{T fx'j) as a function of \z\ for the 
choice of the parameters described in TableUi Dotted, dashed and solid lines correspond 
to aigz = 0, 71 / 4 and tt/2, respectively. 
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Figure 7: BR{t 67) as a function of (f)2 for the choice of the remaining phases described 
in Tahle^ \z\ = tan/? = 10 and Aq = 0. Dotted, dashed and solid lines correspond 

to a.Tgz = 0, 7r/4 and 7r/2, respectively. 
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Figure 8: The extremal values of BR{t e'^)/BR{T — > jj,-^) as a function of \z\ for the 
choice of the parameters described in Table\^ tan/? = 10 and Aq = 0. Dotted, dashed and 
solid lines correspond to aigz = 0, 7r/4 and 7i/2, respectively. 
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Figure 9: Contour plot of BR{fi —>■ ej)/ f{mo, M1/2) minimized with respect to the vari- 
ation o/02 as a function of \z\ and M2/M3 for the choice of the parameters described in 
Ta^/eQ tan/? = 10, = and argz = 0. 
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Figure 10: The CP asymmetries generated in the decays of the sneutrino-inflaton for the 
pattern 0) (left panel) and (0j (right panel) . Details of the scanning procedure are given 
in the text. 
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Figure 11: Minimal values of BR^fi — >• e'y)/ /{rriQ, M1/2) for tan/? = 10 and Aq = 0. The 
dark ring is the region of z in the complex plane for which leptogenesis is successful for 
Ml — 10^^ GeV. The larger ring hounded by thin solid lines show how this region changes 
for Ml = 5 X 10^1 GeV. 
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